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We compute QCD corrections to the production of a tt pair in association with a hard 
photon at the Tevatron and the LHC. This process allows a direct measurement of the top 
quark electromagnetic couplings that, at the moment, are only loosely constrained. We 
include top quark decays, treating them in the narrow width approximation, and retain 
spin correlations of final-state particles. Photon radiation off top quark decay products is 
included in our calculation and yields a significant contribution to the cross-section. We 
study next-to-leading order QCD corrections to the pp — >■ tt^ process at the Tevatron for the 
selection criteria used in a recent measurement by the CDF collaboration. We also discuss 
the impact of QCD corrections to the pp — ^ ^7 process on the measurement of the top quark 
electric charge at the 14 TeV LHC. 
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I. INTRODUCTION 

More than fifteen years after the discovery of the top quark, many of its quantum numbers 
are still not well-measured experimentally. For example, until recently [l| it was possible to have 
a consistent description of "top" quark physics, assuming that the electric charge of the "top" 
quark is Qt = —4/3, in contrast to its usual value Qt = 2/3 [2]. By analyzing tracks of charged 
hadrons to estimate jet charges, the DO collaboration excludes, at less than 2a confidence level, 
the hypothesis that the event sample comes entirely from a heavy quark with the electric charge 
Qt = —4/3. As pointed out in Ref. jsl, a more direct way to measure the top quark charge is to 
study the production of a top quark pair in association with a hard photon. Indeed, to an extent 
that photons are only radiated off the top quarks, the rate for ti^ production is proportional to 
the square of the top quark electric charge. This assumption works well at the LHC, once photon 
radiation off top quark decay products is suppressed, but it fails at the Tevatron because top quark 
pair production there is dominated by qq annihilation. 

The CDF collaboration has recently measured the cross-section of the process pp — )• ti'^, search- 
ing for an excess in events that contain a lepton, a photon, 6-jets and large missing energy J]. 
Using 1.9 fb~^ of data, they observed nine events that they interpreted as due to the pp — )• ti'y 
process. The analysis will be improved by using a larger data sample [5[|, so about fifty ti^ events 
can be expected in the near future. It is therefore important to have a reliable prediction for 
this process, including the possibility of applying selection criteria for final state particles used in 
Ref. 14 1 . 



The authors of Refs. 



1^ 



6| analyzed the potential of the Tevatron and the LHC to study elec- 
troweak couplings of the top quark, focusing on the tiZ and ti'^ final states. If we neglect parity 
non-conservation, the interaction of top quarks with on-shell photons is described by two quantities 
- the electric charge Qt and the anomalous magnetic moment at- Both of these quantities can be 
studied in the ti^ production process. As shown in Refs. P, 0], the best sensitivity to Qt at the 
LHC is obtained if kinematic cuts force the photon to be emitted either in the production stage or 
in the decay stage. Because the non-vanishing anomalous magnetic moment of the top quark cor- 
responds to a dimension-five non-renormalizable operator, it leads to a harder spectrum of photons 
in pp{pp) — )• tt'^. Given enough statistics, it should be possible to study this effect experimentally. 
Regardless of the details pertinent to a particular measurement, once selection criteria are specified, 
the study of top quark electromagnetic couplings becomes a counting experiment which may be 
subject to significant higher order QCD corrections. Therefore, the computation of next-to-leading 
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order (NLO) QCD corrections to tt'j hadroproduction, that correctly incorporates decays of top 
quarks, becomes important. 

NLO QCD corrections to the production of a ti pair and a hard photon in hadron colhsions 
were recently calculated by Duan et al. [7]. This computation was performed in the approximation 
of stable top quarks. While such an approximation gives an idea about the significance of higher 
order QCD effects for the production of ti^, it can not be used to find the magnitude of NLO QCD 
corrections when specific cuts are imposed on top quark decay products. As we explained in the 
previous paragraph, the ability to do this is important for a realistic analysis. In this paper, we 
extend the results of Ref. 3] by computing NLO QCD corrections to pp{pp) — )• + 7, allowing for 
decays of top quarks. We note that radiative decays of top quarks are included into our analysis. 
To calculate one-loop virtual amplitudes, we employ the method of generalized D-dimensional 
unitarity suggested in Ref. ^ and extended to massive particles in Ref. Q]. The current paper 
builds upon the previous studies of tt and tt + j production in hadron collisions, performed by two 



of us 



Id, 



Many technical aspects of the calculation are explained in those references and we 



do not repeat them here. 

When top quarks are treated as truly unstable particles, non-factorizable QCD corrections 
appear [l^. Non-factorizable corrections imply a cross-talk between production and decays of top 
quarks; they can not be described in the narrow width approximation. It is well-understood by now 



12|] that, in many cases, these non-factorizable corrections lead to effects that are suppressed by 
0{asTt/mt), instead of the naive 0{as) expectation for the suppression. Recently, the smallness of 
non-factorizable corrections in reactions with top quarks was confirmed by an explicit computation 
of the NLO QCD corrections to p pipp ) — t- W^W^bh process that included both factorizable and 



non-factorizable contributions 
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IJ]. In what follows we ignore the non-factorizable corrections 



and work in the on-shell approximation for top quarks. 

As a final comment, we note that NLO QCD corrections are known for two other processes 



where the top c 
and on — )• tiZ 



two 



uark pair is produced in association with color- neutral objects - pp 
and to the production of a ti pair in association with one 



20(1 jets, as well as in association with a bb pair j2ll424i|. In all the cases, the NLO QCD 
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corrections are calculated either assuming that all final state particles are stable, or treating QCD 
radiation in top decays incompletely^. Similarly to the ti^ case, removing these omissions may 



^ See however Ref. [2^ where QCD radiation in top quark decays is included in the computation of the top quark 
pair production cross-section in association with one jet by means of a parton shower. 
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become important for precision phenomenology, especially when aggressive cuts are involved to 
separate signals from backgrounds. 

The remainder of the paper is organized as follows. In Section |TI] we describe the setup of the 
calculation and present some results for the case when top quarks are treated as stable particles. 
In Section IIIII we discuss the computation of NLO QCD corrections to radiative decays of top 
quarks. In Section HVl we present phenomenological studies relevant for the Tevatron and the LHC, 
including decays of top quarks. We conclude in Section |Vl Technical details of the calculation are 
described in the Appendices. 



II. PRODUCTION OF A tt PAIR AND A PHOTON: STABLE TOP QUARKS 



We first discuss the case of stable top quarks. To compute the NLO QCD corrections to 
pp{pp) — >■ ttj, we need to calculate one-loop virtual corrections and to account for the emission of 
an additional massless parton. For the calculation of the virtual corrections, we employ the method 



of generalized D-dimensional unitarity suggested in Ref. 



8|]. This method has been used earlier by 



two of us in the computation of hadroproduction of + jet in Ref. {llj. To describe the ti^ final 
state, we can re-use much of that computation. For example, linear combinations of color-ordered 
one- loop amplitudes for ^ tt + 3g [11] give color-ordered amplitudes for — )• tt -|- 2(7 + 7 [26|. 
Similarly, color-ordered amplitudes for —)■ tt -|- gg + 5 [Ul] can be used to construct color-ordered 
amplitudes for — >• tt + qq + ^. All the details of how amplitudes with gluons and quarks are 
transformed into amplitudes with gluons, quarks and a photon are given in Appendix |Al We have 
checked our results for virtual corrections by re-calculating them, for a few phase-space point s by 
using an independent implementation of the Ossola-Pittau-Papadopoulos (OPP) procedure 27l |. 
that we apply to individual Feynman diagrams. The Feynman diagrams are generated with the 
package FeynArts |28]. 



The second, logically distinct part of any one-loop computation is the calculation of real emis 
sion corrections. When integrated over available phase-space, these corrections diverge. Such di- 
vergences must be removed by an appropriate procedure. We use the dipole formalism of Ref. 
extended to deal with QCD radiation off massive particles in Ref. 
calculation can be found in Refs. 1231 . 



29| 



30|. Dipoles relevant for our 



32I ]. In the actual implementation of the subtraction 
procedure, we closely follow Ref. [311] • We have checked that our results do not depend on the 
parameter that restricts the integration over the dipole phase-space; this is a useful way to control 
the consistency of the implementation of the subtraction terms and to improve the efficiency of the 
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FIG. 1: Kinematic distributions in the produetion of a ti pair and a hard photon at the Tevatron, for stable 
top quarks. The bands correspond to the variation of the renormahzation and factorization scales in the 
interval mt/2 < /i < 2mt. Upper (lower) panes show t (7) transverse momentum and rapidity distributions, 
respectively. 



computation [33l |. 

As we mentioned earlier, the calculation of NLO QCD corrections to ttj production in hadronic 
collisions for stable top quarks was reported in Ref . [7.] . When we choose the setup of the calculation 
as close to Ref. as possible, we get good agreement with their results. However, some choices 
made in Ref. 3] ~ for example the use of the electromagnetic coupling at the scale Mz, the use 
of charm and bottom masses in the computation of the partonic channels qg — )• tt'j + g, q = c,b 
and the three degree cut on the opening angle between the photon and the light quark in the final 
state - do not look very appealing to us. For this reason, we decided to present a number of results 
for cross-sections and kinematic distributions which can not be directly compared with the results 
reported in Ref. j?! but which, we believe, correspond to more realistic choices of input parameters 
and better resemble details of experimental analyses. 

We now turn to the discussion of our results for the hadroproduction of a tt pair and a hard 
photon, for stable top quarks. Throughout the paper, we choose the top quark mass nit = 172 GeV 
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35l | for leading and next-to- 



and parton distribution functions CTEQ6L1 and CTEQ6.6M |34l . 
leading computations, respectively. The strong coupling constant is evaluated using one- and 
two-loop running with five massless flavors. To describe emission of the real photon, we use the 
fine structure constant oqed = 1/137. Although this choice should be self-evident because of QED 
Ward identities, we emphasize this fact because in many previous studies of the tt^ production, 
the cross-section was computed with a^j^ = oqedI^z) = 1/128. Using the correct value of the 
fine structure constant is numerically important because it decreases the prediction for the tt'j 
cross-section by about six percent. 

For both the Tevatron and the LHC, we require that the photon is relatively hardpj_^^ > 20 GeV 
and that it is isolated. To ensure that the implementation of photon isolation does not violate 
infra-red and collinear safety, we employ the procedure described in Ref. 
considered isolated and events are rejected unless the condition 

iSpartons \ V 7J / / 

is fulfilled for cones of sizes R that are smaller than R^j = 0.4. In Eq.([T]) R^i is the photon-parton 



36( 1 . The photon is not 



angular distance = ^J (y-y — yiY + (¥'7 — ^iY-, where y^^i (f-y^i) are the laboratory frame rapidi- 
ties (azimuthal angles) of the photon and the parton i, respectively. Also, E±^i is the transverse 
energy of the parton i and E_i_^^ is the transverse energy of the photon. We apply all other selec- 
tion criteria to jets if and only if their separation from a photon exceeds R^ j. A jet reconstructed 

rr 

inside the cone of size R-yj is not subject to selection criteria, see Ref. [36]. As our default, we 
set the renormalization and factorization scales equal to each other and choose them to be equal 
to the mass of the top quark /u = mj. We find the cross-section for pp — )• ttj at the Tevatron 
(Vs = 1-96 TeV) to be 

CTLO = 39.97;j^:^I fb, CTNLO = 37.61°:? fb, (2) 

where the lower value correspond to the scale set to fi = 2mt and the upper value to the scale set to 
/i = mt/2. QCD corrections greatly reduce the uncertainty in the predictions for the cross-section, 
changing it from about thirty percent at leading order to about ten percent at next-to-leading 
order. The cross-section for pp — )■ ttj at the 14 TeV LHC is 

cTLo = 1.96i[jj4 pb, ^TNLO = 2.93l|]:^2 pb. (3) 

The residual scale uncertainty in the NLO QCD cross-section at the LHC is about fifteen percent. 

We note that results for ttj production in Eqs. (l2|3p show significant differences in the QCD 
corrections at the Tevatron and the LHC. At the scale n = nit, the NLO QCD corrections decrease 
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FIG. 2: Kinematic distributions in the production of a tt pair and a hard photon at the LHC, for stable 
top quarks. The bands correspond to the variation of the renormahzation and factorization scales in the 
interval mt/2 < fi < 2mt. Upper (lower) panes show t (7) transverse momentum and rapidity distributions, 
respectively. 

the tt'j cross-section by about six percent at the Tevatron and increase the ttj cross-section by about 
55% at the LHC. It is peculiar that the magnitude of the NLO QCD corrections to ttj production 
at the Tevatron and the LHC is very similar to the magnitude of the NLO QCD corrections to tt 
pair production. While the degree of the correlation between these corrections is perhaps somewhat 
surprising, it can be understood, at least partially, by considering emissions of soft photons which 
must factorize from the production process even when the NLO QCD corrections are included. 
This may also explain why a very similar pattern of QCD corrections was reported in Ref. Q], 
despite the fact that somewhat different input parameters are used in that computation. 

Kinematic distributions of top quark and photon transverse momenta and rapidities at the 
Tevatron and the LHC, are shown in Figs. 11121 No dramatic changes in shapes of these distributions 
are observed. However, the rapidity distribution of top quarks at the Tevatron exhibits some 
interesting features. Indeed, in the process pp — t- ttj, photons can be emitted from both initial and 
final states. Interference of these emissions gives rise to a charge or forward-backward asymmetry 
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of top quarks. We define 

^ ajyt > 0) - ajyt < 0) , . 

* a{yt>0)+a{yt<0)' 

where yt = — In — ] is the rapidity of the top quark in the laboratory frame. As was 

pointed out in Ref. (7|, the forward-backward asymmetry is significant. Calculating it in leading 
and next-to-leading order in perturbative QCD, with the parameters specified at the beginning of 
this Section, we obtain 

At,LO = -17.21°:° %, ^t,NLO = -11.91l3 (5) 

where the central value corresponds to the renormalization and factorization scales set to fi = nit 
and the lower(upper) value to fi = mt/2 and fi = 2mt, respectively. 

It is peculiar that the change in the NLO QCD asymmetry is nearly identical to the size of the 
NLO QCD corrections to forward-backward asymmetry in pp — )■ ttj, computed in Refs. {ill. [isl. fiol . 
In fact, this similarity of NLO QCD corrections to the asymmetries in pp — )• ttj and pp — )• ttj is 
easy to understand, following the observation of Ref. [11] that large NLO QCD correction to the 
asymmetry in pp — )■ ttj is related to the 5% asymmetry in the inclusive rate for pp — )■ tt, first 
computed in Ref. [37]. 



III. QCD CORRECTIONS TO RADIATIVE DECAYS OF TOP QUARKS 

In this Section, we describe the computation of the NLO QCD corrections to the radiative decay 
of the top quark t — )• bWj. Since this is a relatively low-multiplicity process, the calculation of 
virtual corrections is performed using conventional Feynman diagrams, generated with FeynArts 



281 ]. and the Passarino-Veltman reduction [38|. For a few kinematic points, the results are checked 



23], that 



against a computation based upon an independent implementation of the OPP procedure 
is applied to individual Feynman diagrams. 

The real emission corrections to the decay rate are computed using the subtraction formalism 
described in Ref. [s^, in the context of non-radiative top decay t — )• bW. However, the subtraction 
is also applicable to the radiative decay t — )• bWj if one replaces the invariant mass of the VF-boson 
with the invariant mass of the VF-boson and the photon in all the relevant formulas in Ref. [391 ]. 



Specifically, the subtraction term, required to make the real emission contribution t — )• Wb'j + g 
integrable is given by the product of the matrix element squared for the process t — )• Wb'j and the 
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dipole that reads [39|] 

D{pt,pg,pb) = ATiasli^^Cp 



-J-(-^-l-z-ye{l-z)) 

PbPg \l- Z J 



{PtPgf 



The kinematic variables used in Eq.([6]) are defined through 

2 2 

PbPg = ^{'i^ - r) y, ptpg = —{1 - r )(1 - z), 



(6) 



(7) 



where = {pw + p.y)'^ / ml . In Ref. [39|], the dipole in Eq.([6]) is integrated over the unresolved 
phase-space 

1 ymax 

j [dg] D{pt,pg,pb) = M J dz (r^ + z{l - r^)y' J dyy~'{y„,^, - y)-'D{pt,Pg,Pb), (8) 



where e = (4— d)/2 is the parameter of the dimensional regularization, d is the number of space-time 
dimensions and 



{l+rfz{l - z) 



1 + r 



2e 



(9) 



z + r'^{l-z) ' 167r2 r(l-e)Vl-?- 

It is convenient to restrict the subtraction counter-terms to parts of the phase-space that are 



not too far from the singular limits 



33] ■ If this is done, the subtraction terms need to be modified. 



Introducing such a modification in Eq.(6) and integrating over restricted phase-space, we find 
j V^g] D{pt,Pg,Pb) [1 - 6'(1 - a - z)e{y - at/max)] = 



2tt mfT(l - e) 
1 



1 1 /5 



e V2 



2x\ 25 1 / 1 



7 Inr^ 



+ 



2(1 -r2) 



+ 2Li2(l-r^) 



5^ 
~6" 



V 



-21n2a 



'7 ^ a2\ 2(1 -a)r2, 

__4a+— lna+ ^ _ In 



1 — a + r^a 



(10) 



We now present some numerical results for the QCD corrections to the radiative decay of the 
top quark t — )• bWj. We use as{Tnt) = 0.107691, which corresponds to the CTEQ NLO value 
of the strong coupling constant at /i = Mz supplemented with the two-loop running to = mt. 
We take the mass of the VF-boson to be M\\r = 80.419 GeV. We work in the top quark rest 
frame and require the photon energy to be larger than 10 GeV and the opening angle between 
the momentum of the bottom quark and the photon to be such that cos 9b^ < 0.98. With these 
input parameters, we obtain the radiative decay width of the top quark t — )• bW^ at leading and 
next-to-leading order in perturbative QCD 



Flo = 4.48 MeV, Enlo = 3.89 MeV. 



(11) 
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FIG. 3: Kinematic distributions in radiative decays of top quarks t Wbj. Upper panes show normalized 
distributions of the cosine of the angle between the 6-quark and the photon and of the photon energy, 
computed through NLO QCD. Lower panes show ratios of NLO and LO kinematic distributions. 

For the choice of the kinematic cuts described above, the QCD radiative corrections decrease 
radiative decay width by thirteen percent. This is similar to, but somewhat larger than, the 
magnitude of the NLO QCD corrections to the top quark decay width t — )• bW, which decrease 
the decay width t — t- Wb by about eight percent. 

In Fig. [3l we show distributions of the opening angle between the bottom quark and the photon 
and of the photon energy, and ratios of NLO and LO distributions. Shapes of these distributions 
are perfectly described by leading order computations; the NLO QCD corrections provide an overall 
renormalization factor. The distribution of the photon energy shows canonical enhancement of the 
soft photon emission probability at low E^, while the distribution in the opening angle shows a 
collinear enhancement peak at small relative angles between the bottom quark and the photon. 
Suppressing emissions from bottom quarks is important for the analysis of the top quark charge 
that we discuss in the next Section; a simple way to accomplish this is to require that the 6-jet and 
the hard photon are sufficiently separated. 
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IV. PRODUCTION OF A tt PAIR AND A PHOTON: UNSTABLE TOP QUARKS 

A realistic description of the hadroproduction of the tt'j final state requires including decays of 
top quarks. This is important for two reasons: first, it defines realistic acceptances and, second, 
photons can be radiated from the top quark decay products and it may be important to take 
this effect into account. In this Section, we present some results for realistic selection cuts at the 
Tevatron and the LHC. 

The computation described below is performed within the following framework. For both the 
Tevatron and the LHC, we consider the lepton plus jets channel pp, (pp) — )• (t — )• W^{l^i')b){t — >■ 
W~ {jj)b)'y. Top quarks are treated in the narrow width approximation and all spin correlations 
are retained. We include decays of VF-bosons into leptons of definite flavor (e or /i but not both) 
and hadronic decays of VF-bosons into two families of light quarks, that are always treated as 
massless. The VK-bosons are on their mass-shells and we do not consider QCD radiative corrections 
to hadronic decays of PF-bosons. We include photon radiation in the production of a tt pair and 
photon radiation in the decays of top quarks. We note that photons can be radiated by any charged 
particle in the top quark decay process, including the decay products of VF-bosons. 

Before describing the results of the computation, we summarize the input parameters. For 
numerical calculations we use mt = 172 GeV, Mw = 80.419 GeV, the value of the Fermi constant 
Gf = 1.16639 X 10-5 GeV~2 and Fw = 2.14 GeV. We use CTEQ6L1 and CTEQ6.6M 



341, 



35| 



parton distribution functions, for leading and next-to-leading computations, respectively. For 
Os we use one(two)-loop running for leading(next-to-leading) order calculations, neglecting the 
contribution of top quarks to the QCD /3-function. For reference, we give numerical values for the 
top quark width T{t — )• Wb) at leading and next-to-leading order ^] 



= fi_^yfi + 2^) = 1.4653 GeV, Ff = 1.3396 QeV, (12) 



where F^'^^ is calculated with as{mt) = 0.107691. 



A. tt'j production at the Tevatron 

The CDF collaboration has recently measured the production cross-section of the pp — )• ti^ 
process 4]. A new analysis is under way [s] that will extend it to a larger data sample. It is therefore 
of interest to compute the NLO QCD corrections to pp ^ ~^ W^{l^h')b){t ^ ^~(jj)^)7 using 



selection criteria that are employed in the ongoing analysis. Following Refs 



0], we impose 
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FIG. 4: Kinematic distributions of photons in pp {t ^ W''^ {l^ i^)b)(t — > W~ {jj)b)j process at the 
Tevatron. The bands correspond to the variation of the rcnormahzation and factorization scales in the 
interval mt/2 < fi < 2mt. We show the transverse momentum and the rapidity distributions of the photon 
as well as the distribution of the azimuthal angle - rapidity distance R(,^ between the photon and the 
hardest 6-jet. Finally, we show the NLO QCD X-factor as a function of photon transverse momentum 
where factorization and renormalization scales are set to /i = mj. 

constraints on transverse momenta and rapidities of leptons, photons and jets in the process 

p±^i > 20 GeV, p_L,^ > 10 GeV, p±j > 15 GeV; (13) 
\yi\ < 1.1, \y^\ < 1.1, \yj\ < 2. 



In addition we require that there is missing transverse energy in the event E_i 



> 20 GeV 



and that the transverse energy H± is larger than 200 GeV. We define the transverse energy as 
H± = -E^^miss + X]j-^-L,i! where the sum includes the charged lepton, the photon and jets in the 
event. The photon must be isolated from jets R^j = J {y^ — yj)'^ + {ip^ — ipj)'^ > 0.4 and leptons 



R^i = ^ {yj — yiY + {if^ — ipif' > 0.4. The photon isolation is implementec follow 



owing Ref. [36[|; we 



42] with AR = 0.4, 



described it in Section HH We define jets using the A:^-clustering algorithm 141 
and require that at least three jets are reconstructed, two of which are 6-jets. The 6-jets are 
defined as jets that contain 6-quarks from top decays, when partons are clustered according to the 
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FIG. 5: Fraction of events originating from photon radiation in the production pp — tt^ and decay t — > Wbj 
processes, computed at NLO QCD. The renormahzation and factorization scales are set equal to /i = nit- 

jet algorithm. 

With the cuts specified above, we obtain the following results for the total cross-sections at 
leading and next-to-leading order in QCD 

C7LO = 2.85^^:1^ fb, cTNLO = 2.64^0:21 fb. (14) 

The central values correspond to the renormahzation and factorization scales set to fi = mt and 
lower(upper) value to = 2mt{mt/2), respectively. The scale uncertainty in the NLO cross-section 
is reduced by a factor of four, compared to the leading-order one. For fi = nit, the NLO QCD 
corrections reduce the cross-section by about eight percent, similar to stable top quark results 
discussed in Section [III 

In Fig. [3]we show kinematic distributions of photons in the process — t- — )• W^{l^i')b){t — 
{jj)b)'y. We observe a significant reduction in the scale dependence for all kinematic distri- 
butions. The shapes of these distributions do not change much although the photon transverse 
momentum distribution becomes somewhat softer. This is illustrated in the lower-right pane in 
Fig. m where the local i^T- factor K = d(TNLo/do"Lo is shown in dependence of the photon transverse 
momentum. Assuming that the integrated luminosity of 10 inverse femtobarns will, eventually, be 
analyzed by the Tevatron collaborations, tt pairs accompanied by photons with transverse momenta 
as high as ~ 100 GeV should be observable. Li Fig. [5] we separately show fractions of accepted 
events where the photon is radiated either in the production or in the decay process, computed 
through NLO in perturbative QCD. We observe that low-p± photons are produced with comparable 
probabilities in the tt production and decay stages, while photons with high transverse momentum 
P±,7 > 80 GeV are mostly radiated in the production stage. In Figure EJ we show distributions of 
the charged lepton transverse momentum and rapidity, as well as distributions of missing energy 
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FIG. 6: Kinematic distributions in pp {t ^ W'^{l~^i')b){t — W~ {jj)b)'-f process at the Tevatron. The 
bands correspond to the variation of the renormahzation and factorization scales in the interval mf/2 < 
/i < 2mf. We show the transverse momentum and the rapidity distributions of the charged lepton and 
distributions of the missing energy and H±^ . 



and i^T- Distributions of lepton transverse momentum and missing transverse energies become 
softer at next-to-leading order. The reduction in the forward-backward asymmetry for top quarks, 
discussed in the previous Section, is visible in the rapidity distribution of the charged lepton. 



Finally, we remark that the cross-sections shown in Eq. (|14p correspond to our best approxi- 
mation to the setup of the experimental analysis^ described in Ref. [J], and it is interesting to 
estimate the number of ttj events that this cross-section value corresponds to. To do so, we take 
the NLO QCD cross-section shown in Eq. (ll4p . multiply it by a factor of four, to account for the 
possibility to produce e^, final states and multiply by a factor of 0.22 which roughly reflects the 



experimental efficiencies [5|. Final 



y, we multiply by the luminosity 1.9 fb ^ which corresponds to 

4 X 0.22 X 1.9 fb"^ X 2.64 fb w 4. It is 



the data sample analyzed in Ref. Jl and obtain A'^events 

n 

peculiar that in the experimental analysis [4] nine tt^ events were observed. Since a measurement 
that uses a larger data sample is under way, it will be interesting to see what happens to this 



The implementation of photon isolation in the experimental analysis is different from what we use in this paper. 
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difference. Repeating tlie above estimate for a luminosity of 6 ft , we find 14 ttj events. 



B. tt-f production at the LHC 

Next, we study tfie pp ^ (t ^ W^{l'^i')b){t — t- W~'{jj)b)^ process at tfie LHC witli tlie center 
of mass energy of 14 TeV. We apply the following generic cuts that describe detector acceptances 
and the trigger 

p±^^ > 20 GeV, \y^\ < 2.5, R^^b > 0.4, Ryj > 0.4, R^^i > 0.4, 

p±^b > 20 GeV, p±j > 20 GeV, p±^i > 20 GeV, £;±,miss > 20 GeV, 

\yb\ < 2.0, \yj\ < 2.5, \y,\ < 2.5. (15) 



We require that there are two 5-iets anc 
using the fc^-clustering algorithm 
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at least two light jets in the event. Jets are defined 



42] with AR = 0.4. We require 



arge transverse energy 



H± > 200 GeV. The photon isolation is implemented following Ref. 36!]. Using these cuts we 
obtain the cross-sections for pp — t- tt'j production 

^TLo = 74.50™ fb, (7NLO = 138ti fb. (16) 

The central values correspond to the renormalization and factorization scales set to fi = mt and 
lower(upper) value to ^ = 2mt{mt/2), respectively. It follows from Eq. (ll6p that QCD corrections 
are rather large. We saw in Section |lT] that NLO QCD corrections to pp — )• ttj process, evaluated 
in the approximation of stable top quarks, increase the production cross-section by a factor 1.5, 
for = nit. It follows from Eq. ()16p that when acceptance cuts Eg. (USD are applied, the i^- factor 
increases to 1.86 and there is only marginal decrease in the scale dependence. We have checked 
that this increase is related to the radiation of an additional hard jet in pp — )• ttj process and 
that this enhancement disappears when the additional jet is required to be relatively soft. The 
large scale dependence of the NLO prediction for the cross-section is caused by the contribution of 
the quark-gluon partonic annihilation channel that only appears at next-to-leading order. While 
this effect exists at both the Tevatron and the LHC, it gets significant enhancement at the LHC 
due to a much larger gluon luminosity. We note that similar X-factors and large residual scale 
dependence can also be observed in pp — )• tt production at the LHC, when basic kinematic cuts are 
applied to top quark decay products. 

In Figs. I7|8I we show various kinematic distributions at the LHC. Among other things, we 
observe a dominance of the contribution from photon radiation in top quark decays over radiation 
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FIG. 7: Kinematic distributions in pp {t ^ {l^iy)b){t ^ (ji)^)7 process at the 14 TeV LHC, using 
cuts specified in Eq. ljlSp . The bands correspond to the variation of the renormalization and factorization 
scales in the interval rnt/2 < fi < 2mt. We show distributions of the transverse momentum and the rapidity 
of the photon, as well as the distribution of the rapidity-azimuth distance between the photon and the 
hardest 6-jet. We also show the if -factor in dependence on the photon transverse momentum and the 
fraction of events for photon radiation in the tt production and decay stage. 

in the tt pair production, for p±^^ < 60 GeV. As an illustration, we quote results for the NLO 
QCD cross-sections vi^here photon radiation occurs either in the production or in the decay stage 

<,^o^60.9fb, arc"" = 77.2 fb. (17) 

These results correspond to the factorization and the renormalization scales set to the top quark 
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FIG. 8: Kinematic distributions in pp ^ {t ^ W+{l+iy)b){t W~{jj)b)j at the 14 TeV LHC, using cuts 
specified in Ea. (|15p . The bands correspond to the variation of the renormalization and factorization scales 
in the interval mj/2 < /i < 2mj. We show transverse momentum and rapidity distributions of the charged 
lepton as well as distributions of the missing transverse energy -Ex^miss and the transverse energy H±. 

mass, ^ = mt\ their sum gives the total NLO cross-section shown in Eq. (jl6p . We also see that the 
spectrum of emitted photons becomes harder, in contrast to the Tevatron case. The X-factor, as 
a function of the photon transverse momentum is shown in Fig. [71 From Fig. [7] we estimate that 
photons with transverse momenta of up to 350 GeV should be observable at the 14 TeV LHC with 
ten inverse femtobarns of accumulated luminosity. 

As mentioned earlier, the production of ti^ at the LHC can be used to constrain electromagnetic 
couplings of the top quark and, in particular, its electric charge. In what follows we compare some 
results for top quarks with charges Qt = 2/3 and Qt = —4/3. Needless to say that the "top" quark 
with the charge Qt = —4/3 decays through t — )• bW~^ , to respect electric charge conservation. With 
acceptance cuts shown in Eq. (fT5]) . we find 

a^r"^' = 136.8tii; -Slo"'/' = 2431^° fb. (18) 

The central values correspond to the renormalization and factorization scales set to fj, = mt and 
lower (upper) value to = 2mt{mt/2), respectively. These results, together with the cross-sections 
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FIG. 9: Kinematic distributions in pp — ^ ^(7) l^vhhjj + 7 at the 14 TeV LHC for two top quark electric 
charges. The bands correspond to the variation of the renormahzation and factorization scales in the interval 
7714/2 < ^ < 2mt at NLO QCD. 



values shown in Eq. (jl6p . imply that the ratio of Qt = —4/3 and Qt = 2/3 cross-sections is 

-4/3 Qt=-4/3 



Qt=2/3 



1 0/1+0.02 T^NLO _ '^NLO 



a 



Qt=2/3 



1.76: 



"0.01 
-0.02- 



'LO ^NLO ' 

The uncertainty in the ratio comes from the dependence of the production cross-sections for the 
two top quark charges on the renormalization and the factorization scales fi. In Fig. [9] we compare 
basic kinematic distributions for two top quark charge assignments, Qt = 2/3 and Qt = —4/3. 

We note that the dominance of gg annihilation at the LHC seems to suggest that cr{pp — )• tt^) 
cross-sections should scale like the electric charge of the top quark to second power, so that the 
naive expectation for TZ in Eq. p9|) is 7^ = 4. It is obvious from Eq. (|19p that this expectation 
fails. This happens because photons are dominantly radiated by the decay products of t and t (cf. 
Fig. [7]) and this contribution does not scale as Qt- It is also interesting to remark that shifts from 
leading to next-to-leading order and the remaining NLO uncertainties in the cross-section ratios 
shown in Eq. ()19p are quite small, in particular when compared to the corresponding uncertainties 
in the cross-sections, Eqs. (jl6ll8p . In fact, it is easy to imagine that large changes in pp — )• tt^ cross- 
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FIG. 10: Transverse momentum and rapidity distribution of the charged lepton for two values of the top 
quark electric charge in pp tt{'j) -> l^vhhjj + 7 at the 14 TeV LHC, normalized to pp ^ tt ^ l^vhhjj. 
In the latter case the top quark charge is kept at its canonical value Qt = 2/3. 



sections, from leading to next-to-leading order, are not particular to ti^ production and originate, 
rather, from the underlying dynamics of pp — t- ti process. To prove that this assertion is valid, we 
compute ratios of pp — )• tt(7) — )• I'^vhbjj + 7 to — )■ it — )• I'^vhbjj cross-sections, subject to basics 
cuts shown in Eq.([l5]), at the 14 TeV LHC. We obtain 



a 



=2/3 



tt-y 



LO; 



6.33 



-0.26 



X 10-3, NLO, 



Qt=-4/3 
ttj 



hO.3 



LO; 



(20) 



_o.i4 - , — , - ll-2lo:2 X 10-3, NLO. 

It is clear that ratios of cross-sections are significantly more stable against radiative corrections 
and scale variations than the cross-sections themselves. Moreover, these ratios help reduce scale 
uncertainties in kinematic distributions as well. We illustrate this in Fig. [10] where we show 
lepton kinematic distributions in pp — )• ttj for Qt = 2/3 and Qt = —4/3 at next-to-leading order, 
normalized to similar distributions in pp — )• tt. In both cases, basic cuts shown in Eq. (|15p are 
applied. It is striking that for Qt = 2/3, the ratio datt^/dau is essentially constant for a large 
range of kinematic parameters, while for Qt = —4/3 the relevant spectra appear to be harder. 
However, in both cases the scale uncertainty of the ratio is much smaller than the scale uncertainty 
when pp — )■ ^7 and pp — t- tt are considered separately. It is clear that, in addition to the scale 
uncertainty, other uncertainties such as in as and in parton distribution functions cancel to a large 
extent in the cross-section ratio [431 ], making it an interesting observable to study at the LHC. 



Although the ratio of Qt = —4/3 and Qt = 2/3 cross-sections shown in Eq. (jl9p appears already 
large enough to distinguish between the two electric charge assignments, one can make it even 
larger. Indeed, as we already mentioned, at the LHC top quarks are mostly produced in gluon 
collisions. Hence, in the production stage of the process, photons are radiated mostly by top 
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FIG. 11; Fractions of pp ttlj) — ^ l^vhhjj + 7 events at the 14 TeV LHC with photon radiated in the 
production and in the decay, in case when RDS cuts are apphed. The renormahzations scale is set to ^ — mt. 

quarks. If we manage to reduce the contribution from the radiation in the decay, we will have 
an observable that is sensitive to the electromagnetic couplings of the top quark. To reduce the 
probability that the photon is radiated in the top quark decay, we impose the following cuts: 

• we determine the 6-jet jb that forms the smallest invariant mass with the charged lepton 
and require large transverse mass of that 6-jet, lepton , photon and the missing energy 
'n^±{jb^l'y P±miss) > 180 GcV. The transverse mass here is defined as 

m]_{jblT,P±,miss) = (^yjp'iijbh) + m^Ubh) +P±,miss^ - {p±{jbh) +P±,miss)^; (21) 

• the remaining 6-jet is combined with the two hardest light jets; it is required that the invariant 
mass of these three jets is close to the top quark mass 160 GeV < m{jjjb) < 180 GeV; 

• to suppress photon radiation from leptonic decays of VF-bosons, we require m±{£j] -E^miss) > 
90 GeV; 

• to suppress photon radiation from hadronic decays of VF-bosons, we require that there are 
two light jets in the event whose invariant mass is close to the mass of the VF-boson, 70 GeV < 
m{jj) < 90 GeV. 

We will refer to those cuts as the "radiation-in-the-decay-suppression" (RDS) cuts and we 
emphasize that they are applied in addition to cuts shown in Eq. (|15|) . Applying RDS cuts, we find 
that the leading order tt^ production cross-sections reduces by about a factor of three compared 
to the case when only generic cuts Eq. ljlSp are applied 

CTLO = 23.39j:^:| fb, CTNLO = 26.7li| fb. (22) 
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Compared to cross-sections shown in Eq. (jl6p . there are significant changes in the NLO result as 
well since the X-factor becomes much smaller when RDS cuts are applied. Kinematic distributions 
are shown in Figs. Ill|12[ The degree of suppression for photons originating from top quark decays 
can be seen in Fig. [11] where the distribution of the photon transverse momentum and the angular 
distance between the photon and the hardest 6-jet R'yj^ are displayed. We observe that with the 
RDS cuts more than ninety percent of the total cross-section is due to photons radiated in the 
production of a tt pair with about five percent coming from top quark decays. 

Finally, we apply the RDS cuts to compute the cross-section for the production of the top quark 
with the exotic charge Qt = —4/3. We find 

a^r'^' = 72.62tiri fb, a^lZ-'^' = 76.9^ fb. (23) 

where the central value corresponds to ^ = nit and the lower (upper) value to /x = 2mt and 
/i = mi/2, respectively. We determine the ratio of the cross-sections for the two charge assignments 
and find that it increases 

Qt=-4/3 Qt=-4/3 

T^LO _ ^LO _ 9 1 n+0-05 ^NLO _ '^NLO _ o qq+O.OS fc)A\ 

'■^RDS — Qt=2/3 ~ "^■-'■^-0.04' '■^RDS " Qt=2/3 ~ ^■°°-0.12' l^^J 
'''lo "^NLO 

compared to Eq. (fT9]) where only basic cuts are applied. 

We can now estimate if it is worth applying the RDS cuts. We denote by C the luminosity 
required to separate Qt = —4/3 from Qt = 2/3 at the 3a level with the cuts in Eq.([T5]) and by 
>Crds the same quantity when the RDS cuts are applied in addition. The two quantities are related 
by the following equation^ 

_c_ ^ (7^RDs - 1)^ 

C^DS CtQ*=2/3 (7^ - 1)2 ^ 

We can use Eqs. (|16|18|19|22|23|24|) to compute the ratio of the required luminosities at leading 
and next-to-leading order in perturbative QCD. Interestingly, because the X-factors for the two 
types of cuts are so different, we find that the required ratios of luminosities differ by a significant 
amount 

£ I 1.98 ±0.02, LO; 
J-RDS I 1.12 ±0.08, NLO. 

It follows from Eq. (l26p that once next-to- leading order effects are accounted for, the application of 
RDS suppression cuts becomes much less important since a factor of two gain in luminosity gets 



^ We only consider statistical errors. 
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FIG. 12: Kinematic distributions in pp tt{j) — )■ l^vhhjj + 7 at the 14 TeV LHC, for two electromagnetic 
charge assignments of top quarks with RDS cuts. The two lower panes show normalized distributions, 
to emphasize differences in shapes. The bands correspond to the variation of the rcnormalization and 
factorization scales in the interval mt/2 < /i < 2™^. 

reduced to 0(10%) gain. We also find that kinematic distributions are not very sensitive to the 
top quark charge; for illustrative purposes we show some distributions that have some sensitivity 
to Qt in Fig. [T2j The largest effect is present in the lepton transverse momentum distribution that 
becomes harder when the top quark charge increases. Also, the rapidity distribution of the charged 
lepton becomes more central, but this effect is not very significant. 

V. CONCLUSIONS 

In this paper, we describe the calculation of the NLO QCD corrections to the production of a 
ti pair in association with a hard photon at the Tevatron and the LHC. This process is of interest 
for direct studies of the electromagnetic couplings of the top quark such as its electric charge and 
its anomalous magnetic moment. In a recent measurement, the CDF collaboration detected nine 
ti'y events in pp collisions using an integrated luminosity of 1.9 fb^^. Our best estimate, which 
includes NLO QCD corrections, realistic acceptances and photon radiation in the production and 
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decay stages of a tt pair, is that 4.4 it 0.2 events should have been observed with the integrated 
luminosity of 1.9 fb^^ and an efficiency of 22%. It will be interesting to compare our results with 
an analysis of a larger data sample, that is under way Q]. 

For any hadron collider process, an accurate prediction requires at least next-to-leading order 
QCD computation and this is what we set out to do in this paper for the ttj final state. For 
stable top quarks such a computation was performed in Ref. [7i]. However, for practical studies 
of ti^ production, the result of Ref. [t^ is not sufficient since, in general, large fraction of isolated 
photons comes from radiative decays of top quarks. It is important to incorporate top quark decays, 
including radiative ones, into a unified framework that also includes higher-order QCD corrections, 
and this is what we accomplished in this paper. 

We studied the pp{pp) — t- tt{'~f) — t- l^vhbjj + 7 process both at the Tevatron and at the 14 TeV 
LHC, including effects of the NLO QCD corrections. We found that, in general, the QCD correc- 
tions are small at the Tevatron and are large at the LHC. However, we also observed that these 
QCD corrections are very similar to the NLO QCD corrections to pp{pp) — )• tt processes suggesting 
that ratios of these cross-sections datt^/dau can be theoretically predicted with higher accuracy 
than the two cross-sections separately. 

We found that about fifty percent of all photons in the tt'j events are radiated off the top quark 
decay products, both at the Tevatron and the LHC. This fraction increases with the decrease in 
the photon transverse momentum reaching approximately eighty percent for p_\_^'y ~ 10 GeV at the 
LHC. Since photon emission off the top quark decay products is a background to measuring elec- 
tromagnetic couplings of the top quarks, it is important to apply selection criteria that suppresses 
these contributions. Designing cuts to suppress photon radiation in top quark decays and applying 
them to pp —7- tt^ process at the LHC, we observe a significant change in the i^- factor compared 
to the i^-factor computed for basic cuts. This feature emphasizes the importance of fiexible im- 
plementation of the radiative corrections to processes with unstable particles, where kinematics of 
the decay products must be accessible. 

To have a concrete model where the electric charge of the top quark is different from its Standard 
Model value, we have studied the case when Qt = —4/3 and the "top quark" decays into the 6- 
quark and the VF~-boson . We have used this charge assignment to investigate the possibility 
of measuring Qt by studying ttj production at the LHC. If the analysis is performed at LO QCD, 
we find that designing cuts to suppress the QCD radiation from top quark decays benefits the 
analysis. However, when the same analysis is performed at NLO QCD, suppressing QCD radiation 
off the decay products of top quarks becomes less important because, when basic cuts are applied, 
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the tt'j production process receives large corrections at next-to-leading order in perturbative QCD. 

Finally, we note that the production of photons in association with a tt pair can be studied at 
the 7 TeV LHC. For the selection criteria as in Eq. (jl5[) . we find the leading and next-to- leading 
order cross-section for pp — )• tt(7) — )• l^ubbjj -1- 7 to be 15 fb and 26 fb, respectively. If 5 fb""*^ are 
indeed collected at the 7 TeV LHC by the end of the year 2012, we estimate that about 500 events 
with high energy isolated photons, large missing energy, an isolated lepton and two 6-jets should 
be observed at the LHC. Among these events, there will be a few truly spectacular ones, with 
the tt pair accompanied by a very energetic photon. We look forward to studies of the pp — )• tt'j 
process at the LHC in the coming years. 
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Appendix A: Gauge-invariant decomposition of scattering amplitudes with a photon 

In this Appendix, we present the gauge invariant decomposition of helicity amplitudes that we 
used in the calculation reported in this paper. Our starting point is the computation of the ttj 



ladroproduction process reported in Ref. 



ll[ |. From that reference, we know primitive amplitudes 
2a] for partonic processes such as — )• ttggg, — )■ ttqqg etc. and we would like to turn them into 
amplitudes that describe production of a tt pair, quarks, gluons and a single photon. We do so by 
constructing linear combinations of tt + gluons + quarks primitive amplitudes in such a way that 
non-Abelian contributions cancel. Because the photon can also be radiated in the decay of the 
top quark, amplitudes for tt pair production without photon radiation are also required. Those 
amplitudes can be found e.g. in Ref. [10], but we present them here for completeness as a special 
case of the tt^ amplitudes. 

We begin with the amplitudes for leading order processes for tt^ production. At leading or- 
der, two partonic initial states gg and qq contribute to the cross section. We write the color 
decomposition of the corresponding matrix elements in the following form 

M'^'^gg ^ tt7) = g^ V2eQt ^ {T^-^T'^^^l -^''''i^t, ^, 2^, (as),, (^4),), (Al) 

0-GS2 
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In Eq. (|M|) we use 



jii^is jiijia 

«4 «2 _/\r «2 «4 



Q 



2i , 3,-, 4,; 5^) = 5„ 2t, 3,-, 4,) + ^ ^^-^(If, 2*, 3,-, 5^, 4, 



(A2) 



(A3) 



We note that tree partial amplitudes A}''^'^'^ in Eqs. (IAllA3p involve only quarks and gluons, i.e. 
no photon. The factor ^/2 in Eq. (|Al|A2p compensates for the similar factor in the color-stripped 



Feynman rules 



441 ] that we use to compute amplitudes with quarks and gluons. Because we require 



photon, rather than gluon, emission amplitudes, we must remove this factor. 

To obtain amplitudes for the tt final state, without photon radiation, we use the following set 
of rules: 

• remove 5^ from partial and primitive amplitudes; 

• set \pieQt — )• 1 and Qg, Qf ^ 0; 

• set auxiliary parameter k — )■ 0. 



These rules can be applied to Eqs. (|Al|A2IA3P as well as to all other photon-emission amplitudes 
that we present in this Section. 

We now describe amplitudes required for the next-to-leading order QCD computation. For the 
real emission corrections we require four partonic channels gg — )■ ttjg, qq — t- ti^g, qg — >• tt'jq and 
qg — )• tt'jq. The last three channels are related by crossing symmetry. For this reason, we present 
the color decomposition for gg — )• ttg^ and qq — )• tt^g. We find 

M''^\gg^ti^g) = gl V2eQt J] (T'^-3r'^^4T'^^6)| 5„ 2^, (^3),, (^4),, (ag),), (A4) 



M'^^'^'iqq ^ tt^g) = glV2eQt 



+ (r"«)1^5i:Br(lf,2i,3,-,4„6,;5,) 
+ ]^(^''^)S<5&^r(lf,2*,35,4„6,;5,) 



1 



+ ^(^"^)S'^l^r(lf,2*,35,4„6,;5,: 



(A5) 



where 



2i, 3,-, 4„ 63; 5^) = A'^'^'ilt, 5g, 2t, 3,-, 4„ 6^) + ^ A'^'^'^ili, 2* , 3,-, 5„ 4„ 6^), (A6) 



2i, 3,-, 4„ 63; 5^) = 5„ 2t, 63, 3,-, 4,) + ^ A'^'^'i^, 2t,%, 3,-, 5„ 4,), (A7) 



Qt 

Qg 

Qt 
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^s^'^ih^ 2i, 3q, 4g, 6g; 5^ 



+ 1 / ^tree(i_^ g^^ 3*, 3,-, 4,) + A'^'^^U, 6„ 5„ 2* , 3,-, 4, J 
+|^^*-'=(lf,63,2i,3,-,5g,4,)), 
Sr^(li,2i,3,-,4„63;5^) = 5„ 2,, 3,-, 6^, 4,) + ^ ( ^^-''(If, 2^, 3,-, 5„ 6„ 4,) 



\-ti -yi -tj -yj 7 ^y/ ' 

+^*'-^<^(li,2i,35,63,55,4g)y 



(A8) 



(A9) 



As long as we are interested in tt^ amplitudes, k should be set to one. For tt amplitudes without 
photon emission, k, should be set to zero. 

The matrix elements of the virtual corrections are constructed from primitive amplitudes. For 
the partonic channel gg — )• tt'j we find the color decomposition 



M-^^'igg ^ ttj) =gtV2eQtY. E 



(2^a.3ya.,)n^virt(i_^2i,(a3),,(a4)g;5,) 



+Tr(r»''3r«<'4) ^IJ^vf 2t, (^3)a, (cT4)a; 5.) 



(AlO) 



with 



^i,b (It; 2j, 3g, 4g; 



^2,a (1*; 2j, 3g, 4g; 



^2,b*(lt> 3g,, 4g; S-yy 



N,A'^'^'\li,5g,2t,3g,Ag) - ^|^(^^'W(lf,4„5„3g,20 



5g, 4g, 3g, 2t) + 4^, 3^, 5^, 2*) j , (All) 

(If, 5„ 2t, 3„ 4g) + 1^ (1,-, 2t, 3g, 4g, 5g) 



+4'^'/'] (If, 2^,3,, 5g,4g) + (If, 2i, 5„ 3„ 4,) 

+4'[^/'](lf, 5^,2^,3,, 4,))|, 

^L'W (If, 5„ 2t, 3gAa) + (^"^'f'' (If, 5^, 4„ 3^, 2t 

4g, 5g, 3g, 2t) + ^L'W(lf, 4c;, 3^, 5^, 2 



+^ ( ^"^'t'kli, Sg, 3g, 2^,4^) + ^L'W(lf, 3^, 5^, 2t, 4^, ; , 



■^i,b*(lt; 2t, 3g, 4^; 5^;, 



(A12) 



(A13) 



(A14) 



where ^'^'l^l and are regular primitive amplitudes as defined in [26]. The sum over / in 



Eq. ()A12p and Eq. (jA14p includes all quark flavors. The color decomposition of the matrix element 
for the process qq — t- ttj follows Ref. [45,] and is given by 



M''"\qq^tt^)=gtV2eQt ^ 5|<5|^r;r(lf, 2^, 3^, 4,; 5 



i=a,b,c,d 



+ 5i^^j^^j;f(lf,2,,35,4,;5, 



(A15) 
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FIG. 13: Definition of primitive amplitudes witli fom" quarks. Figure (a) defines Aa whiclr includes all 
topologies where both external fermion lines enter the loop. Figures (b) and (c) define At and Ac where 
either the top quark or the light quark line enter the loop, respectively. Finally, Figure (d) defines amplitude 
Ad that includes topologies with a closed fermion loop. 



with 



^™ai^ti 2t, 3q, Aq] 5^ j 



^l,b*(lt) 2t, 3q, Aq-, 5-) 



^l,c (It) 2t, 3q, iq] 5-, 



^l|d*(lt' ^i?' ^q'l ^1 



^^a(lf, 5g, 2t, 4q,, 3g) — ^^a(lt, 2t, 4g, 5g, 3q)^ , 

2k + 1 



A(lf, 5g, 4,, 3,, 2t) + A(lf, 4g, 5g, 3g, 2t 



1 



A(lf, 5^, 2t, 3q, Aq) - A(lf, 2t, 5g, 3g, 4g) 

+Ac{lt, 2t, 35, 4g, 5g) + Acil-t, 5g, 2t, 3g, 4g)|^ , 

J2 (^iih, ^, 2t,3q, Aq) + 2i, 3,-, Aq)) . 



f 



(A16) 



4g, 3,-, 5^, 2t) - ^MU, 4g, 5^, 35, 2t)) , (A17) 



(A18) 
(A19) 



The primitive amphtudes ^a,b,c,d correspond to different topologies depending on which fermion 
lines enter the loop, see Fig. 13. We note that the contribution of the amphtude vanishes 
after its interference with the tree amplitude is computed; for this reason, we do not present it 
here. 



Appendix B: Formula for the production cross-section 



In this Appendix, we present the formula that we use to describe radiative corrections to top 
quark pair production in association with a photon in the narrow width approximation. In the 
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narrow width approximation, the differential production cross-section is given by 

dc7 = dau-,dBt,xdBt,Y + da^f [dBt,X',dBt,Y + dBt,xdBt,Y-y\ , (Bl) 

where Bt^x{-i) is the branching fraction for either radiative t — t- X + 7 or non-radiative t ^ X top 
quark decay. The above equation is valid to all orders in the strong coupling constant. We expand 
through 0{a\) and denote 

j-nLO Jr<5NLO p<5NL0 

rLO _ '^^ i^X R.5NLO _ Q-*^ i~^X ^, _ 1 o^t /-RO^ 

- -pLO ' - pLO ' X - - ^ pLO ' ^^^^ 

where i G is the total decay width of the top quark at leading order and T^^^ is the 

QCD correction to the total top quark decay width. We obtain 



da = X ( + + x d4? 

+ daf ( AB\^^ m\% + AB\'^ AB\%:) 

+ d4o( dB^L^o dB^^ + dBl^ dBff.^ + dBff^ dBl% + dS^;0^ de^^Lo) 

Since we also treat the M^-boson in the narrow width approximation, photon radiation in the top 
quark decay can be further decomposed into photon radiation in top quark (t — > Wb) and W 
{W ff') decays. The decays of the VF-bosons are treated at leading order in perturbative QCD. 
We therefore write 

d^U. = ^ + ^ ^g^, ^ e [LO, 6NL0], X = bff. (B4) 
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